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Abstract. We explore observables in a lattice Universe described by a recently found 
solution to Einstein field equations. This solution models a regular lattice of evenly 
£■ — distributed objects of equal masses. This inhomogeneous solution is pcrturbative, 

exact at second order in a small parameter, and expands at a rate exactly equal to the 
(J one expected in a dust dominated Friedmann-Lemaitre-Robertson- Walker (FLRW) 

model with the equivalent, smoothed, energy density. Therefore, the kinematics of 
5— i both cosmologies are identical. Looking at the behaviour of the redshift and angular 

i distance, we find a condition on the compactness of the objects at the centre of each 

cell under which corrections to the FLRW observables remain small. Nevertheless, 
we show that, if this condition is violated, namely, if the objects are too compact, 
our perturbative scheme breaks down as far as the calculations of observables are 
t-H concerned, even though the kinematics of the lattice remains exactly identical to its 

FLRW counter-part. This may be an indication of an actual fitting problem, i.e. a 
situation in which the FLRW model obtained from lightcone observables does not 
correspond to the FLRW model obtained by smoothing the spatial distribution of 
matter. 
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PACS numbers: 04.20.-q, 04.20.-Cv, 98.80.- Jk 
1. Introduction 



It has long been recognised that calculating observables in an inhomogeneous Universe 
could be quite challenging. Indeed, our usual description of the geometry of the Uni- 
verse on large scales relies on Friedmann-Lemaitre-Robertson- Walker (FLRW) models 
in which the distribution of matter is assumed perfectly homogeneous and isotropic. 
As a result, light rays propagate in a homogeneous medium filled with matter and are 
therefore sensitive to the Ricci curvature of spacetime. In the real Universe, on the other 
hand, at least in the late stages of its evolution, matter is clumped into virialised objects 
with large (almost) empty regions between them, and light therefore travels mainly in 
empty space, where its behaviour is dominated by the Weyl curvature of spacetime, 
rather than the Ricci curvature. This raises the natural question to whether the FLRW 
approximation is suitable to calculate observables in the late time Universe: under 
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which conditions is it possible to replace Weyl curvature along the line of sight with an 
'equivalent' Ricci curvature? And what is this 'equivalence' all about? The problem is 
therefore what has been dubbed a 'fitting problem' [1], that is: how do we replace the 
real Universe by an FLRW idealisation? In principle, this leads to a backreaction issue: 
the idealised FLRW given by observations will not, in general be the FLRW model for 
which the actual matter density has simply been smoothed in space. In the standard 
model, this is usually accounted for via the Dyer-Roeder equation [2] , but recent works 
tend to show that there might be better ways to model the effect [3, 4]. This issue has 
been addressed by a certain number of authors lately; see e.g., [4, 5, 3, 6]. Usually, they 
use a 'realistic' model of structure formation, either through cosmological perturbation 
theory or N-body simulation, to evaluate the impact of inhomogeneities on the propaga- 
tion of light; others have used Swiss-Cheese models [7, 8, 9, 10], mainly to see whether 
inhomogeneities along the line of sight could account for the Dark Energy phenomenon. 
Our take on the problem is slightly different. In particular, we are not trying to address 
the Dark Energy problem; we simply try to understand under which conditions, in the 
controlled environment of an (almost) exact solution to Einstein field equations, the be- 
haviour of null geodesies can safely be approximated by the null geodesies of an FLRW 
model. 

In a previous paper [11], we proposed a lattice solution to Einstein field equations 
made of equal masses M separated by a comoving distance L. This solution is exact 
at order M/L, and can be expanding or contracting. We proved that this solution is 
kinematically equivalent to an FLRW model with a dust matter content having an en- 
ergy density equal to the one obtained by smoothing the lattice distribution, i.e. M/L 3 . 
This was interpreted as supporting the usual fluid approximation in cosmology, and was 
consistent with results obtained by others [12]. Such a solution, being kinematically 
equivalent to an FLRW model, is the ideal setting to study observables: if observables 
only deviate slightly from the observables in the kinematically equivalent FLRW Uni- 
verse, then we might say that this FLRW model is a good fitting model of our Universe. 
On the other hand, if observables exhibit large variations compared to the one in the 
analogue FLRW, that might be a sign that the fitting problem in cosmology has to be 
studied extremely carefully, since this shows an explicit example of a spacetime that 
does not lead to the same, averaged FLRW model when considered kinematically, or 
along its null geodesies. 

The propagation of light in lattice models has been studied before [13, 14] in the context 
of Lindquist- Wheeler models [15]. These models are only approximate, and the propa- 
gation of light through the boundaries between cells is not fully controlled, leading to 
important differences in the results, depending on the approximation scheme used [14]. 
Our model, on the contrary, is only approximate in power of M/L, and not on the way 
cells are glued together (much like the model developed in [16]). Therefore, propagating 
light through the lattice is not a problem. The only limitation will come from the fact 
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that we can only trust the solution for a small range of redshifts. By solving Sachs 
equations at order M/L, we will see that at this order the equations for the shear and 
the isotropic expansion decouple. This implies that at order M/L, a priori, the Weyl 
curvature of spacetime does not play any role in the calculation of the distance/redshift 
relation. Then, the formal solution appears to be equivalent to the solution in the ana- 
logue FLRW model, plus small corrections of order M/L. Nevertheless, by studying 
carefully the order of magnitude of these corrections, we will prove that they remain 
small only under specific conditions on the compacity of the objects forming the lattice: 
when these objects are too compact, the perturbative expansion breaks down, and the 
differences between observables in the lattice and in the analogue FLRW model cannot 
be simply evaluated by using perturbative methods. That might indicate a tension be- 
tween the fitting models constructed kinematically and observationally. 

The paper is organised as follows. Section 2 will briefly present the solution described 
in details in [11]. Section 3 details the calculations of the redshift and angular distance 
in the lattice Universe. For the sake of clarity, Section 3 focusses only on the analytical 
expressions, while their physical interpretation, numerical calculations, and discussion 
of the effect of inhomogeneities on the propagation of light, especially in regards of the 
fitting problem and of the Weyl focussing, are left for Section 4. Finally, section 5 will 
be a conclusion. 



Throughout the paper, the signature of the spacetime metric will be (— , +, +, +), and, 
unless otherwise specified, we will work in units of c = G = 1. 



2. The lattice solution 



In a previous paper [11] we found that an infinite, regular, cubic lattice of masses 
M spaced by a comoving length L, in an otherwise empty Universe without any 
cosmo logical constant, is described by the following metric in a synchronous comoving 
coordinate system: 



9oo = ~ 1, 9oi = 0, 



1 + 2e 



GM 




8tt ct 2GM 
3 L + Lc 2 



/(*) + 



2nc 2 t 2 
3L 2 



- ^c¥^/(x), 



Lc 2 tJJ ^ 

up to order GM/Lc 2 , and where / is given in Eq. (3) below. This perturbative, 
approximate solution, is justified phenomenologically by the fact that a lattice of galaxies 
separated by the typical intergalactic distance today have GM/Lc 2 ~ 10~ 8 . In the above 
metric, e — ±1 corresponds respectively to an expanding or contracting lattice. In the 
following we shall restrict ourselves to the expanding lattice only (e — 1). This metric 
solves Einstein equations without cosmological constant up to order M/L, with the 
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source term associated to the lattice, namely: 

neZ3 v y v y 

where the point-masses are described by a three-dimensional Dirac comb. In this case, 
the function /, accounting for the anisotropics of the gravitational field created by the 
masses, reads: 

ft \ = \^r = cos ( 2 fnx) cos ( 2 fpy) cos ( 2 fqz) 

/W 7T ^ |n| 2 "tt ^ n 2 +p 2 + 9 2 " ' [ ) 

where N 3 stands for N 3 \ {(0,0,0)} and similarly for Z 3 . This solution presents 
singularities at the position of the masses as a natural consequence of modelling them 
with a Dirac distribution. However the solution can be regularized in the UV, e.g. 
by describing instead the lattice of masses by a three-dimensional lattice of peaked 
Gaussians of width r\. In this case the non-vanishing part of the stress energy tensor 
now reads 

Too = S(x)S(y)S(z) : (4) 

with: 

x 1 (x-nL) 2 

and similarly for S(y) and S(z). This source term can be written in Fourier space as: 

rp M ^ £ in . x _*W 

neZ 3 

Then, the above metric, Eq. (1), is still a solution to Einstein equations provided the 
function / is now replaced by the following regularized function f v : 



1 ^ e L2 j. 



2n in.x 



^8 e g cos { T nx) cos ( T pyJ cos (-^gsj , , 

7i ^ n 2 + p 2 + q 2 

An important feature of both the source terms presented in Eqs. (2) and (6) is the fact 
that the infinite wavelength mode (the zero mode n = of the Fourier expansion of the 
source), behaves as an homogeneous comoving density M/L 3 . This Friedmannian-like 
component for the source precisely produces in Eq. (1) the FLRW metric of a flat, dust- 



dominated Universe, when expanded in powers of Ho = y 87r/3 x M/L, and up to order 
Hq, see [11]. On the other hand, finite wavelength modes of the source distribution 
have no Friedmannian counterparts and account for its inhomogeneity. These modes 
generate the terms proportional to f v and its derivatives in Eq. (1). The metric of the 
lattice Universe is thus given by the FLRW metric, truncated at the correct order, plus 
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corrections coming from the inhomogeneity and anisotropy of the source distribution. 
As a consequence, the observables in the lattice Universe will partially match the ones in 
FLRW, up to the degree of approximation considered here, but with corrections coming 
from the anisotropic terms in the metric. Therefore, our main task in this paper will be 
to assess the amplitude of these corrections with respect to the FLRW behaviour at the 
level of observables quantities, such as redshift and angular distance. 

The validity of the approximations made depend essentially on two dimensionless param- 
eters, M/L and rj/L; their ratio is simply the compactness of the objects. As mentioned 
above, the numerical values of M and L can be chosen to be of the order of magnitude 
of the mass of a galaxy and of a typical intergalactic distance: 

M ~ 10 n M o and L ~ 1 Mpc, 

where M Q is a solar mass. The choice of rj, the spread of the objects, is more diffi- 
cult to make. A natural cut-off is provided by the Schwarzschild radius of the masses: 
rj ~ 2M. However, in the cosmological context, it is more relevant, as well as much 
easier for numerical purposes, to choose a factor r] that matches with the size of the 
objects considered here, i.e. the size of a typical galaxy, which is much larger than 
its Schwarzschild radius. In that case, a good choice would be: r\ ~ 10 kpc, so that 

_ TT 2 (n 2 +p 2 +q 2 )t 1 2 

rj/L ~ 0.01. With this value of rj, the UV-regulator e makes the various 

sums converge quite quickly, and we shall thus restrict ourselves to a large but finite 
number of terms (typically of order 200 3 ) in the expression of the function f v , when we 
will deal with numerical studies in the next section. 



Note that this choice of numerical values is only indicative. In the rest of the pa- 
per, we will keep M/L fixed, but r\/L will be allowed to change, since the main result 
of the paper is the fact that observables only remain close to their FLRW analogue 
provided: 

M «o ( i)xm 4 , ( 8) 



L v ' \L, 

i.e. provided the objects are not too compact. If this bound is not satisfied, the 
perturbative expansion used in this paper fails to be well-defined, and other, non- 
perturbative techniques should be used to address the fitting problem. Note that the 
choice of 77 mentioned earlier, rj/L ~ 0.01, that corresponds to the value expected for a 
galaxy-like object, is marginally violating condition (8), since it corresponds to the case 
M/L ~ {ji/Vf. 



In the following we shall also use some shortcuts for the metric, which we write for- 

Mr , M 

T og ah + -j- 



mally as a power series g a b = f] a b + \l ^f$9ab + y~S 2 g a b, where rj a b is the Minkowski 



metric, and 5g a b and 5 2 g a b correspond respectively to the order a/ M/L and M/L parts 
of the metric in Eq. (1). Accordingly, the Christoffel symbols are decomposed into 



pa / M rpa _j_ M r2pa 

be \l L be ' L be' 
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3. Observables: analytical expressions for redshift and distance 

In order to evaluate observables, we investigate the properties of null geodesies in the 
lattice Universe. In the following, k a refers to the dimensionless four velocity of a null 
ray k a = dx a /d\, where A is an affine parameter along the ray and has the dimension of 
a length. In Section 3.1 we solve perturbatively the geodesic equations in order to find 
the redshift z(X) as a function of the arrival direction of the photon at the observer, see 
Eq. (22) below. In Section 3.2, we then use and solve perturbatively the Sachs optical 
equations for a bundle of light rays in the lattice Universe. This enables us to compute 
the expansion 6 of the light-bundle, and thus the angular distance as a function of the 
affine parameter r^A), see Eq. (40) below. In this section, we will provide analytic 
expressions for z(X) and r^(A). The discussion of the amplitude of each term in these 
expressions is left for Section 4. We shall work first in natural units and re-establish 
dimensionful constants c and G at the end of the calculation. 



3.1. Geodesic equation for light rays 

The geodesic equation reads: 
Dk a dk a 



together with k a k a = 0. We look for a pertubative solution of the form: 

[M M /m 3 / 2 \ 
r (A) = v a (\) + yj -e(A) + T C(\) + O (^J . (10) 

The zeroth order is simply v a = 0, where a dot refers to d/dX. We pick up the solution 
v a = (l,v x ,v y ,v z ) with v\ + v y + v 1 = 1, so that the i/'s shall represent the arrival 
direction of the photon at the observer (see below). The first order equations (in a/ M/L) 
then read: £ a + 5T^ c v b v c = 0. From the metric Eq. (1), this explicitly reads: 

These equations are solve by: 




provided the boundary conditions are set such that the observer is located at A = 
and to = 0. Then A is increasing with the cosmic time and is negative along the past 
light cone of the observer. We will take similar conventions for ( a , so that we will have 
& a (Ao) = (l,v x ,v y ,v z ), and the v^s thus give the direction of the light ray arriving at 
the observer, as announced. This solution can then be used to deduce the expression 
for the look-back time to the leading order: 
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Notice that the expression for the look-back time does not involve the usual present 
time to- Actually the fact that to = comes from the constraint k 2 = at first order, as 
one might check directly. Accordingly, the present time to was set to zero by convention 
in the metric Eq. (1). The reason behind the convenience of this choice is that there 
is no notion of absolute time in the lattice Universe, as the dynamical evolution of the 
lattice solution into the past cannot be followed until it reaches any singularity. Indeed, 
denoting L phys (t) the physical distance between the masses (basically equal to L times 
an effective scale factor, see [11]), the approximation M / L phys (t) <C 1 would break down 
at a certain point in the past, because L phys gets smaller and smaller. 

Assuming that the observer is located at x % = 0, we also have: 

x l (\) = v i X + 0{y/M/L). (15) 
The second order equation then reads: 

( a + 25T a bc v b C + S 2 T a bc v b v c = 0, (16) 

where it is enough, at this order, to replace t by A and x l by v l \. Then, the temporal 
component of the geodesic equation reads: 

C°(A) + A^(x(A)>V = (17) 

where summation is meant on repeated Latin indices, and f v was given by Eq. (7). 
One may derive a similar equation for ( l and also check that the constraint k 2 = is 
respected to this order. As far as the redshift is concerned however, we only need k° 
and thus (°. Noticing then that: 

~ [/„(x(A)) - W,(x(A)K] = A3y/,(x(A))vV, (18) 



the differential equation for £ u is solved by: 

\„ 

C°(A ) - C°(A) = [/,(x(A)) - A^x^y]* + 



14ttA 2 



3L 2 



(19) 



A 



With the boundary conditions explained previously, namely Ao = and C°(0) = 0, we 
finally get: 

C°(A) = ^ + [/,(x(A)) - A^(x(A)K] Q A , (20) 
The redshift can now be deduced using: 

where the subscripts S and O refer respectively to the source and the observer, and 
where we used that fundamental observers, such as the masses themselves in our case, 
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have u a = (1, 0, 0, 0). We also used that k°(X ) = 1 as a consequence of the normalisation 
chosen. We thus get the following law for the redshift: 



GM 8n X 



This expression coincides with its FLRW counterpart for a flat, dust-filled Universe 
without any cosmological constant, up to order Hq oc GM/L 3 c 2 . The anisotropies only 
show up in the term into brackets, and at second order in ^JM/L. This formula is 
further discussed in Section 4.1. 



3.2. Sachs optical equations for a bundle of light rays 

A bundle of light rays is described in General Relativity by the Sachs optical equations 
[17]: 

JQ 1 

+ e * + &a = S R = --R ah k a k b (23) 
aX I 

d ^ + 2a6 = S c = C abcd m a k b m c k d , (24) 

where A is an afline parameter along the light ray, 6 and a are, respectively, the isotropic 
expansion and complex shear scalars of the ray bundle. We have set the vorticity u = 
because this is always a solution of the equations for an astrophysical isotropic source. 
Alternatively the Sachs equations can be written for the angular diameter distance 
given by: 



Then they read: 

1 d 2 r A 
T\ dX 



2 



+ aa = S R (X) (26) 



£ + = < 27 > 

The vector field m a in Eq. (24) is a complex null vector field (we denote by m a its complex 
conjugate) such that m a m a = 0, m a m a = 1 and m a k a = 0. It can be decomposed into its 
real and imaginary parts by introducing two spacelike unit vectors n\ and n^, spanning 
the screen space orthogonal to the instantaneous direction of propagation of the light 
rays and parallely transported along the ray: m a = (n\ — in?,) / V2. For the basis in the 
screen space, we can choose, in (t,x,y,z) coordinates, as long as v z ^ If: 

»t = U -P=j, -VT^l) (28) 

| Similar expressions can be used if v z — 1. The system being symmetric by exchange of the axes of 
symmetry of the lattice (aligned with the coordinate axes), it is enough to consider one particular case; 
the other ones can be straightforwardly deduced from it. 
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which correspond to the unit vectors of the spherical basis, the unit vector v = (v x ,v y ,v z ) 
being the radial one. A calculation using the metric Eq. (1) then yields explicit 
expressions for the Ricci and Weyl type sources Sr and Sc- 

S R (\) = -\R A n l = ™ (~ + A/„(x(A))) (30) 

ScW = CWft'**** = -*f l«'*-fU*W\ (31) 

L 1-vj 

where u = v A V and w = u A v, and where we used to this order t = A and 
x l = v 1 X. We note that the sources for the expansion and the shear, Sr and Sc, 
are of order M/L. As a consequence, aa in Eq. (26) is formally of order M 2 /L 2 , so 
that the shear and the expansion do not couple to each other up to the degree of 
approximation considered here. In order to see this more explicitly, let us look for a 
perturbative solution to Eqs. (26, 27, 30, 31), and expand r A = + \pf;r^ + 

and a = + \J~^a^ + j^cr^- Then the zeroth order obeys Sachs equations in an 
empty Universe: 

" r f + *«V0) = o (32) 



r (°) d 2 \ 

1 A 



ivr + 7^^t^ ' (33) 

' A 

whose solution is = and = —A, as in Minkowski spacetime, for the appropriate 
integration constants, and for spherical sources (for which the initial shear is set to zero). 
First order equations are then found to be: 



<Pr<P 



(34) 



d 2 X 
daW a W 

The first order equation for ta gives = A\X + B\ and thus only renormalise . In 
the following, we will therefore write r^(A) = —A + 77?"^ (A). The second equation is 
also solved by = 0, provided the sources are considered to be spherical. Finally, the 
second order equations read: 

1 d 2 r { ? 4tt 



A d?\ L 2 



A/,(x(A)) (36) 
^ (2) , 9 ° {2) K- W j 2 /,(x(A)) 



Hence, as announced, the shear and the expansion (or the angular distance) indeed de- 
couple. More generally, from the form of Sachs equations and the fact that their source 
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is of order M/L, it is clear that the shear and the expansion (or angular distance) would 
couple only at the fourth order in y/~M~/ L. An analysis of this coupling would require 
the knowledge of the source terms, and thus of the metric of the lattice Universe, up to 
M 2 / L 2 included. This solution is unfortunately out of reach, however, because of the 
non-linearity of Einstein field equations, that translates in our case into a non-trivial 
mixing of Fourier modes; see [11]. 

The fact that the shear and the expansion decouple means that for a Universe where 
M/L is small, the Weyl curvature does not have a priori any noticeable effect on the 
distance-redshift relations. This is highly unexpected, as one would have thought that, 
the lattice Universe being mainly empty, the propagation of light would be mostly in- 
fluenced by the Weyl curvature. In fact, the decoupling is a formal consequence of 
the expansion in powers of a/ M/L, and thus holds only if this expansion remains a 
valid approximation. We will show in Section 4 that when the lattice Universe gets 
too inhomogeneous (i.e. when, for a fixed M/L, we consider a more and more peaked 
distribution for the masses, that is, smaller and smaller rj/L), the expansion in powers 
of sjM/L actually breaks down, in the sense that second order quantities in M / L 
happen to be of order 0(1). In this case it is not true any more that shear and angular 
distance decouple, and the Weyl focussing gets back into the game, as it is expected in 
a highly inhomogeneous Universe with large and almost empty regions. 

Before going into this discussion, however, we first derive the analytical expression for 
the angular distance. Eq. (37) is solved by: 

rf{X) = £d\'£ X" - A/„(x(A"))) d\". 

The integration is best done using the expression for f n found in the first line of Eq. (7). 
This yields 

, 2n GM A 3 
r A (X) = - X + 



3 Lc 2 L 2 
GM \ / //. • ;;An.v\ 



ne 



where v = (v x ,v y ,v z ). This expression only holds if n.v ^ 0. Let us denote S v the set 
of triplets (n,p, q) G Z 3 such that n.v = 0. The full result is then: 

, 2irGM A 3 

GM ^ (iL + vrAn.v 
Lc 2 V vr 2 (n.v) 3 

GM 2ttA 3 ^mV . . 

neSv 



e L 
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This expression is real valued, and can be further simplified to get, in the end: 

2tt GM A 3 



r A (\)= -A + 



3 Lc 2 L 2 



1+ J2 e 

(n,p,q)eT> v 



+ 



2GM WjV 



7T Lc 2 
Z=4 



E 



X 



E 



=i 



neN3\© 
COS 



-A- 



+ 



Lsin(^f^) 



(40) 



(V.U;) 2 7T (V.U/) 3 

Here the uj's are the following triplets: 

Ui = (n, p, g) , u 2 = (n, -p, -q), u 3 = (n, p, -g) , u 4 = (n, -p, g) , 

and V v = {(n,p,q) E N 3 : 3/ G {1, 2, 3, 4}/u ; .v = 0}. This means that the first sum 
is over all the triplets that cancel one at least of the u^.v, whereas the second sum is 
over all the other triplets. We note that the first two terms in the angular distance in 
the first line of Eq. (40) coincide with their FLRW counterparts for a flat, dust-filled 
Universe without any cosmological constant, up to order A 3 , or equivalently up to order 
Hq oc GM/L 3 c 2 . The other terms however shall be characterized more precisely in 
Section 4.2. 



4. Observables: properties, numerical results and discussion 

So far we derived analytical expressions for the redshift and the angular distance in the 
lattice Universe, Eqs. (22) and (40). Now, we would like to study in more details the 
'small' differences between these observables in the lattice, and the same observables in 
the FLRW model with the same averaged energy density. 

4-1. Redshift 

The expression for the redshift Eq. (22) contains Friedmannian terms coming from the 
zero mode of the source, given by: 



and anisotropic terms given by: 

anisotropic^ = ™ x [f^)) - A9 i /„(x(A))u i ] J , (42) 

coming from the contribution of the other modes. We note that these expressions are 

GM 14ttA 2 I GM /~8n X 



only valid for a limited range of A such that j^^jj- < y y 3 Z- Numerically, this 
gives 

!£< 0.2^1, (43, 
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or also, in terms of redshift (putting back into Eq. (41)), z <C 0.6. This is the natural 
limitation of our perturbative approach. 

Naively, the effect of anisotropies on the redshift must be small for small value of |A| 
(i.e. small redshifts), since in this case the Friedmannian term in ^ M/L dominates 
over the smaller terms or order M/L, provided that M/L <C 1. At larger redshifts (or 
larger |A|), the anisotropies must be again subdominant since they behave at most as 
M/L x A (the function f v being bounded), while the Friedmannian-like term goes as 
M/L x A 2 . However, these conclusions might be affected in some appropriate interme- 
diate range for A, if the absolute value of f v and dif^v 1 is large. Although this is a 
non-trivial question due to the complicated form of f v and its derivative, it is already 
clear from the expression of f v that its magnitude depends on rj/L, i.e. on the width of 
the sources rj compared with the separation L between the sources. In other words, the 
compactness of the lattice plays a major role in fixing the amplitude of the corrections 
to the pure FLRW behaviour. It is therefore of importance to evaluate how and its 
derivative scale with the compactness rj/L. 

The simplest term in Eq. (42) reads 

ir 2 (n 2 + P 2 +q 2 )r, 2 

/,(0,0,0) = - 6 2,2, 2 ■ ( 44 ) 

7r , n 2 + jtr + q l 

As a matter of fact, it is difficult to track analytically its magnitude as a function of 
rj. Numerically, however, we have found the empirical law: ^(0,0,0) ~ 0.73(r]/L)~ im . 
Hence, we see that this correction stays small if 0.73(r]/L)~ 105 <C 1 ^ . Rounding 
the figures, and using Eq. (43), this condition reads rj 3> 1.2M, which always holds 
since the masses need to be at least as large as their own Schwarzschild radius. There 
are two other terms in z amsotropic . The first reads f v (x(\)), and basically behaves as 
/^(0,0, 0) in order of magnitude, since we can consider that the product of the three 
cosines involved are on average of order unity. The last term reads Xdif v v l , and thus 
behaves as (considering again cos ~ 1 and sin ~ 1) 

_ir 2 (n 2 +p 2 +q 2 ) v 2 

»,//~4 X 0(l)x E <'^+^f>;7^ ^. (45) 
L n 2 + p 2 + q 2 

For a given arrival direction of the photon, the |i> l |'s are random but of order unity or 
smaller. Thus we shall bound Xdif v v l by the following sum 

_ir 2 (n 2 +p 2 +q 2 ) v 2 

\xt i^,r X {n + p + q)e ^ 

*dif v v l < 16- x \ ^ ' . 46 

L n 2 + p 2 + q 2 

Numerically, this sum behaves as ~ 0.1(L/r]) 2 ; see Fig. 1. As a consequence, the 
corrections due to inhomogeneities remain small with respect to the M/L Friedmannian 
term provided that 

M «0. 3 m 4 , (47) 



L \L 
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Figure 1. Numerical result: Log-Log plot of the sum found in Eq. (46), 

/ (n + p + q J, as a function of 



S = E(„,p,g)£N2 ^( nV * + ^ + ^) e_ 

77/L, and its fit (continuous curve). In the range of interest for r]/L, it behaves as the 
power law given in the text. 



all calculations done, and using Eq. (43). For the parameters we chose, namely 
M/L = 10 -8 and r]/L = 10~ 2 , this condition is not fulfilled and we should expect 
large deviations to the FLRW redshift. In fact it is not so (see Fig. 2), because the 
bound we derived here is only indicative, in the sense that Xdif v v l is typically much less 
than its superior bound given in Eq. (46). Still, this indicative result roughly determine 
a subspace in the parameters for which deviations to FLRW cosmology are large§. 



Interestingly however, one might also show, using the above results, that the correc- 
tions due to inhomogeneities (of order M/L) remain small with respect to the leading 
order effect in a/ M/L if 

l 3 (!) • (48) 

which is quite close to the previous bound. This shows that large corrections (in the 
redshift) to the FLRW behaviour are typically associated with a breakdown of our per- 
turbative approach, since then order M/L quantities become comparable to \J M/L 
ones. A very similar conclusion will also arise from the analysis of the behaviour of the 
angular distance as a function of ry; see next subsection. 



These bounds Eqs. (47, 48) (and see also Eq. (54) for the angular distance) between 
the three length scales of the problem, namely M, rj and L, mean the following. If 
these conditions are met, meaning essentially that the masses are not too much peaked, 

§ On evaluating numerically the exact expression for the redshift, we have found that z anlsotro P' lc ^ 
z flrw wncn M/L ~ 1000 (rj/L) 4 . This correspond to the condition for the whole anisotropic 
contribution to be of order the FLRW one. 
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then the corrections to the observables || with respect to the FLRW observables remain 
small, and this shall be taken, in this case, as an argument in favour of the use of the 
cosmological principle in standard cosmology (bearing in mind, of course, that we only 
discuss an approximate solution to GR, with a high degree of symmetry). When these 
conditions are not met, however, it does not mean per se that the FLRW solution is not 
valid. It only says that the corrections to the fitting FLRW Universe are so large that 
the very foundation of our approach to the lattice Universe, namely a series expansion 
in powers of y/M/L, cannot be trusted any more. Still, it gives a quantitative estimate 
of when non-linearities must be taken into account in an inhomogeneous Universe in 
order to get sensible results. Moreover it shows that the compactness of the sources 
is a critical parameter impacting the behaviour of the observables in inhomogeneous 
cosmology, regardless of how small the lattice parameter M/L might be. One must 
insist, again, on the fact that the lattice is kinematically identical, on average, to the 
FLRW model with the same energy density. In that sense, there is no backreaction 'a 
la Buchert' [18, 19], irrespective of the compactness of the objects. Nevertheless, our 
results tend to indicate that there might be 'lightcone backreaction': the fitting FLRW 
reconstructed from observables might be different from the one reconstructed by simply 
averaging out spatial inhomogeneities. Only a non-perturbative treatment could answer 
this question. 

All the previous analysis for the redshift has been checked against direct numerical 
calculation of 2 a ™ so * ro ^ c , displayed below. In Fig. 2 we plot ( z Lattice - z flrw^ ^flrw = 
anisotropic i z flrw ag & f unction f t h e FLRW redshift for M/L = 10" 8 , rj/L = 0.01, 
and for some random v l, s. The corrections are less than one percent, although we have 
checked that the discrepancy with respect to the homogeneous Universe increases with 
the degree of inhomogeneity (i.e. for decreasing r//L), as expected in the light of the 
previous discussion. The amplitude of the effect does not depend significantly on the v l, s 
selected. In order to perform numerical calculations, we have truncated the sum defining 
the function f v or its derivative, and summed over non-null triplets (n,p, q) G {0, N} 3 
with iV = 200, which is enough to ensure convergence with such a value for r)/L. The 
noisy character comes from the fact that the photon reaching the observer has to escape 
the local gravitational well from where it is emitted. Notice that such a noise induces a 
natural spreading of the luminosity distance cIl(z) around its main Friedmannian value, 
which is not imputable to the peculiar velocities of the sources, but is rather an intrinsic 
effect coming from the existence of local gravitational wells, and which moreover couples 
to the expansion, i.e. it is a kind of ISW effect. For this reason, such a spreading of 
the Hubble-Sandage diagram shall actually be shared by any inhomogeneous models. 
The amplitude and the amount of spread for the redshift can be evaluated in our case, 
although it stays below observable limits for realistic values of the parameters. 

|| At least for the redshift and the distances. 
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Percentage change in redshift between lattice Universe and FLRW 
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Figure 2. Numerical result: Plotof<5z/z = (z 



y FLRWy z FLRW ag & function f 



y FLRW 



in pcrccnts. These quantities are defined by Eqs. (41, 42). Here M/L = 10 



■q/L = 10 2 , and only the first 200 3 terms of the sum are considered. 



4-2. Angular distance 

The equation for the angular distance, Eq. (40), shows two types of corrections to the 
FLRW behaviour, depending on whether the triplets (n, p, q) G N 3 involved are poles 
(such that u^.v = for some I) or not. Let us first focus on the first correction to 
FLRW, reading: 



-A + 



2tt GM A 3 
~3~ T&l? 



E 



+ 



(49) 



where we recall that V v = {{n,p,q) G N 3 : 3 / G {1, 2, 3, 4}/u/.v = 0}, see Eq. (40) 
and below. This set is given by the intersection of the plane whose normal is v, with 
the lattice of masses. The sum over V v is thus very difficult to perform in general. 
However, we can expect the sum to be negligible for general (random) values of the 
v %, 8, since then we expect T> v to be sparsely distributed in N 3 . The sum involves many 
terms only for some very specific values of v corresponding to the symmetries of the 
cube. Hence we can evaluate that the sum is always less than its particular value for a 



photon propagating along the axes of the lattice, for example, v x 



0, v z = 1, in 



which case V v = {(n,p,Q) : (n,p) G N^}. Then the sum can be computed analytically 
in terms of the Jacobi elliptic function: 

1 / / \ 2 

4 



1 + 



E ■ 



ir 2 (n 2 +p 2 )r) 2 
772 



1 



v e 



- v 2 /l< 



(50) 



where u(t) = 1 + 2]T™t n2 . Numerically, this behaves as ~ 0(1) x (r]/L)- 2 . This 
means that, at most, and in very rare cases, the corrections to the A 3 term can get quite 
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large, in which case this significantly reduces the range of validity of the perturbative 
scheme followed here: 

2nGM\X\ 3 /n (LV FT 
l A l>^T^^V x0 " 7177' ( 51 ) 



3 Lc 2 L 2 \r]J 11 \LJ V M' 



instead of the natural range ^ <C 0.2 J found in Eq. (43). For an inhomogeneous 



M 

Universe with rj = 10~ 2 L, this reduces the range of attainable redshifts to already unin- 
teresting values of order 10~ 2 , thus illustrating again the inadequacy of the perturbative 
expansion when evaluating observables in the lattice Universe. However, this calcula- 
tion only holds for very specific values of v. We do not expect such large corrections in 
general. Still, this term reduces in any case the range of applicability of the perturbation 
scheme followed here. 

The second correction to FLRW in Eq. (40) is more interesting. This expression in- 
volves pole-like terms proportional to l/(u;.v) 2 . Although u^.v cannot vanish strictly 
by definition of V v , it might still get very large in general. When such a 'quasi-pole' 
l/(u/.v) 2 3> 1 occurs, it basically dominates over all the other terms in the sum in 
Eq. (40), in which case the angular distance simplifies to: 



, 2ir GM A 3 



1+ £ a" 

(n,p,q)eV v 



2 GM e 

7T LC 2 (V.UJ 2 

where only the largest pole given by some triplet (n*,p*,q*) and a specific value for is 
considered here. In deriving the previous expression, we also used that cos ^ 27rA ^ ,u '* j ^ \ 

and sin ^ 27rA ^/ u '* j ^ 27rA ^' u '* since v.u^ <ti 1. We note that the angular distance is thus 
typically smaller than the one in the kinematically fitting FLRW Universe since A < 0. 
This was confirmed by direct numerical calculation including all the terms in the sum; 
see below. It can even happen that the angular distance becomes negative, if the -i/'s are 
such that there exists triplets for which the pole is large enough (i.e., essentially when 
the light ray arriving at the observer came very close to one of the masses of the lattice 
along its path). This can only happen when the second order term above (in M/L) gets 
as large as the zeroth order term —A. Such an unphysical result must therefore not be 
seen as a catastrophe per se, but rather as a signal that our perturbative scheme breaks 
down a at certain point. 

In order to get a quantitative estimate of the range of validity of our perturbative 
expansion, we have thus studied the probability distribution of the magnitude of the 
'quasi-pole' terms l/(u;.v) 2 as a function of i]/L. Our method was the following. 
For given random values for the i> l 's on the unit sphere, we looked for the maximal 
value of l/(u ; .v) 2 for all / = 1, . . . , 4, and all triplets (n,p,q) in {0, ...,N} 3 such that 
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(n,p,q) 7^ (0,0,0). We repeated the calculation for a large number (1000) of random 
values for the i> l 's, and thus ended up with a large enough distribution of the maximal 
value of the 'quasi-pole' as a function of N. We thus took the median of this set to 
find the typical value of the largest 'quasi-pole' as a function of N. Now varying N, we 
found numerically that the typical size of the largest 'quasi-pole' as a function of N goes 
like ~ 22 N 4A . The main correction Sr A /r A to Friedmann's law for the angular distance 
then reads, using the previous result and Eq. (52): 

5r A 44 GM 41 _W_ 

— ~ zr-z-N e l 1 . (53) 

7T Lc 2 

In order for the perturbative expansion to remain a valid approximation, we thus require 

< 1. Rounding the exponent to 4, the function N e ~& is maximal for 

N = y/2L/(nr]) and its maximal value reads (r)/L)~ 4 /12.8. Rounding again the figures, 
we thus find that our perturbative approach to the calculation of observables in the 
lattice Universe only holds if: 

(54) 

all calculations done. This bound defines a restricted range of validity of our perturba- 
tive approach, independent of A, besides the natural limitation Eq. (43) of the Taylor 
expansion in powers of a/ M/L. If the free parameters satisfy Eq. (54) above, then it is 
very likely that the pole-like terms do not blow up. In this case then, the corrections 
to the FLRW angular distance will be small. On the other hand, if (54) is not satisfied, 
it is very likely that a large 'quasi-pole' shows up, in which case the corrections are 
large and can even lead to unphysical results such as a negative angular distance. The 
relevance of the above bound has been checked against full numerical computation of 
the angular distance. It turns out numerically that the equation ^ ~ 10 (£) is slightly 
overestimated, and that the transition between the two limiting behaviours between 
small and large corrections to FLRW occurs more precisely around ^ ~ 0.1 (£) ; see 
Figs. 3 and 4 for an illustration of the two cases of small and large corrections. 



Interestingly, this roughly corresponds to the bound found for the redshift^f, in par- 
ticular concerning the exponent, 4. More generally then, we can conclude our analysis 
of Section 4.1 and 4.2, by saying that the general equation 

^«0(l)x(|)\ (55) 

relating the lattice parameter M/L to the compactness of the sources rj/L is the critical 
condition for the corrections to FLRW to stay small, while, under this condition, the 
perturbative expansion followed here stays under control. Although we only deal with 
an approximate solution for an idealized lattice, this result strongly suggests that the 

% Although, as explained in Section 4.1, the strength of the corrections to Z FLRW have been 
overestimated in the analytical treatment. 
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effects of inhomogeneities upon observable quantities in more realistic Universes are to 
stay small provided that the above bound is met. 

What happens in a more inhomogeneous Universe where Eq. (55) is not satisfied is 
a more difficult question. We saw that it implies a priori large corrections to FLRW. 
However, in this case, it turns out that the perturbative expansion breaks down, leaving 
the possibility for higher order terms to compensate for second order corrections, and 
more generally, to compensate between each others. This goes beyond the scope of 
this paper, since this study would require a solution to Einstein equations up to higher 
orders, and a solution to the Sachs equations that couple shear and isotropic expansion. 

We saw that when the Universe is too inhomogeneous (in the sense of the above con- 
dition), the perturbative expansion blows up. We would like to show, in the next 
subsection that this can be related to the fact that the Weyl focussing cannot be ne- 
glected any more, even at order M/L, when the objects become too compact. In a way, 
this will establish the relation (54) as the criterion to determine whether or not, in a 
given Universe, Weyl curvature plays any role in the propagation of light. 



Angular distance as a function of lambda in lattice and FLRW Percentage change in angular distance between lattice and FLRW as a function of redshift 

1000r 1 1 i 0.3 i 1 1 1 1 1 




Figure 3. Left: The angular distance as a function of the affine parameter A < in 
the lattice Universe (dashed line) and FLRW (continuous curve), for M/L = 1CT 8 and 
r)/L = 0.03, i.e. M/L « 0.012(?y/L) 4 , and random v"s. The order of magnitude of 
the effect does not depend on the v l, s selected. The plot is based on 1500 calculated 
points, and the cutoff for the sums is N — 200. The two lines are indistinguishable. 
Indeed, the relative difference is less than 0.4%; see right panel. 



4-3. Effect of the Weyl curvature; break-down of the perturbative expansion 

The shear is given by a = where is given by the differential equation Eq. (37). 

The equation is formally solved by: 

a(2)(A) = _ W A2 K-^ 2 /,(x(A)) dA (56) 
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Angular distance as a function of lambda in lattice and FLRW Percentage change in angular distance between lattice and FLRW as a function of redshift 




Figure 4. Left: The angular distance as a function of the affine parameter A < in 
the lattice Universe (dashed line) and FLRW (continuous curve), for M/L = 1CP 8 and 
i)/L = 0.017, i.e. M/L w 0.12(?7/L) 4 , and random v l, s. The order of magnitude of the 
effect docs not depend on the v l, s selected. The plot is based on 500 calculated points, 
and the cutoff for the sums is N = 200. The angular distance in the lattice Universe 
is significantly less than the FLRW distance, due to the negative contribution of the 
'quasi-poles', as discussed in the text. The relative difference is around 6%; see right 
panel. This illustrates the sharp transition in the behaviour of the observables as a 
function of M/L x (L/t?) 4 . 



Provided n.v ^ 0, a calculation shows that: 

i GML 2 x - e —^e—tr- 
G = ~ n 2 L(l - v 2 ) Lc 2 A 2 ^ n 2 (n.v) 3 X 

A A 2 \ 

1 + 2i7r-n.v + 2tt 2 — (n.v) 2 J (iU z - W z f . (57) 

where U = v A n and W = U A v. The case n.v = leads to terms proportional to 
A. This expression is only defined when v z ^ 1; if v z — 1, a similar expression could 
be found by changing the basis in the screen space (cf footnote before Eqs. (28) and 
(29)). The exact expression does not matter much here. What is important is that we 
thus find similar terms in the shear and in the angular distance, especially the poles and 
'quasi-poles'in n.v. We saw in the previous section that when such 'quasi-poles' blow 
up, that is, when they become numerically of order L/M 3> 1, which is very likely if 
the bound ^ <^ 0(1) X (£) is not met, then the corrections to the angular distance 
due to inhomogeneities become very large and the perturbative expansion breaks down. 
It is thus very interesting to note that the shear itself blow up around the same time 
and that it is no longer of order M/L. In this case, the term aa appearing in the Sachs 
equation for the angular distance is not of order M 2 /L 2 any more, and therefore cannot 
be neglected. Thus, the bound Eq. (55) is also critical as to assess whether the shear 
and angular distance decouple or not at first orders. This shows that, as expected, 
the Weyl curvature plays an important role in very inhomogeneous Universes violating 
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this bound. As a matter of fact the precise point at which the shear blows up cannot 
simply be read out of Eq. (57), because this expression was obtained by neglecting the 
contribution of the angular distance to the equation for the shear; but we saw that this 
angular distance blows up when the bound Eq. (55) is not satisfied, therefore, a full 
solution for the shear must take into account an of order one. This emphasises the 
fact that, in order to have a good understanding of what happens to observables in the 
lattice when the bound Eq. (55) is not met, one must solve the full system of coupled 
Sachs equations, without any perturbative scheme, even though the metric can still be 
written perturbatively. This is beyond the scope of the current work and is left for a 
future study. 

5. Conclusion 

In this paper, we showed that in a lattice Universe kinematically equivalent to an FLRW 
model with the same averaged energy density, observables would not deviate significantly 
from the ones in the FLRW model provided the compactness rj/L of the object was 
related to the mass of the objects M/L via: 



This relation can thus be understood as a criterion to decide whether or not the fitting 
problem is a problem at all: if this relation is satisfied, observables are almost the 
ones on the kinematically fitting FLRW model. Otherwise, if objects are too compact, 
perturbative estimates are not reliable, and it is impossible to say how observables relate 
to the fitting model without a fully non-perturbative treatment of the propagation of 
light. This should not be understood as a claim that there is a genuine fitting problem 
in our 'real' Universe, but simply as a warning that there exist spacetime configurations 
such that, despite the fact that the solution of Einstein field equations remains close to an 
FLRW configuration, observables might significantly deviate from their FLRW analogue; 
see [20] for a similar result in FLRW plus perturbations. Indeed, we managed to link this 
break-down of the perturbative expansion to the fact that the Weyl curvature behaves 
non-perturbatively when the bound (58) is not satisfied, leading to big corrections to 
the angular distance at second order in a/M/L, via the effect of the shear. It would be 
interesting to check if similar conclusions can be drawn in more 'realistic' configuration 
of the matter distribution. In particular, provided one could find a satisfactory way 
of estimation the compactness of cosmological objects, it would be interesting to see 
whether or not the bound (58) is satisfied in our Universe, using galaxy surveys and/or 
N-body simulations. A future work must also present the numerical solution of the full 
system of Sachs equations non perturbatively, in order to avoid the limitations of the 
bound (58); only such a solution will allow one to decide what happens when this bound 
is not satisfied. The results of [21] imply, through statistical arguments, that there is no 
fitting problem in cosmology, even for a matter distribution made of point masses. Our 
model, with a high degree of symmetry, seems to indicate otherwise. For example, the 




(58) 
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'quasi-poles' seem to act only in one way, decreasing the angular distance compared to 
its FLRW counter-part. Nevertheless, it is actually impossible to conclude and one will 
need the full, non perturbative, solution in order to do some statistics and to tests the 
results of [21] with the lattice solution presented here. 
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